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ABSTRACT 

Constant wall temperature /homogeneity in wall temperature is the need of various lab-on-chip 
devices employed in biological and chemical investigations. In this work, a novel and simple 
way of maintaining constant wall temperature is proposed. A diverging microchannel along with 
conjugate effects are utilized towards this end. Both measurements and three dimensional 
numerical simulations are undertaken to prove the design. The investigation has been carried out 
over a large parameter range (divergence angle: 1-8°; length: 10-30 mm; depth: 86-200 |um; 
solid-to-fluid thickness ratio: 1.5-4.0, and solid-to-fluid thermal conductivity ratio: 27-646) and 
input conditions (mass flow rate: 2.5-5.5 ml/min, heat flux: 2.4-9.6 W/cm 2 ) which helped in 
establishing the finding. It is observed that a nearly constant wall temperature condition can be 
achieved over a large parameter range. A model is also proposed to arrive at the design 
parameter values. The method is demonstrated for Polymerase chain reaction (PCR) where we 
are able to maintain temperature each of the three stations within +1 °C. The finding is therefore 
significant and can be employed in both single and multi-stage processes requiring different 
constant wall temperatures with a fine resolution. 



1. INTRODUCTION 

The primary objective of this study is to propose a method for maintaining constant wall 
temperature in microdevices. In several biological applications, it is required to maintain the 
cells above/ below the room temperature implying the need for heating/ cooling the microdevice. 
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For instance, Hung et al. [1] designed a microfluidic cell culture array for studying the human 
carcinoma cells in which a constant temperature of 37 °C was required to be maintained. Wan et 
al. [2] studied proliferation and migration of tumor cells under a controlled temperature of 37 °C. 
Polymerase chain reaction (PCR) is a widely employed technique to amplify a particular DNA 
sequence and is used in a number of applications like pathogen detection [3,4] and hereditary 
disorder diagnosis [5]. The entire process takes place in three steps: denaturation, annealing and 
extension for which the different stations are required to be maintained at constant temperatures 
of 95, 55 and 72 °C, respectively. Selva et al. [6,7] required a certain temperature gradient in 
order to control the movement of bubbles and droplets by means of thermocapillary effects. 
These examples underline the need for developing an effective and efficient system which could 
maintain the microchannel walls at a desired temperature. 

When the temperature of the sample is to be raised, as in the above examples, one can either 
supply a uniform heat flux at the boundaries, or maintain the boundaries at a higher temperature. 
As is well known, for hydrodynamically and thermally fully-developed flow in simple 
geometries, a uniform heat flux leads to a linear increase in surface temperature, and a constant 
difference between the bulk fluid and wall temperatures. The presence of a lateral temperature 
gradient in the fluid however implies that the local fluid temperature can exceed the safe 
operating temperature limit. Maintaining the boundaries at a constant temperature is therefore the 
safer alternative where the local fluid temperature will necessarily be less than the wall 
temperature (in the wall heated case). At the conventional scale, the constant wall temperature 
condition is typically maintained by employing phase change. Maintaining constant temperature 
is however not simple, more so when extended to microscale. Moreover, only certain 
temperatures (and not any desired temperature) can be achieved with the phase change technique 
due to the limitations in the material that can be employed for this purpose. 

Hung et al. [1] and Wan et al. [2] employed an incubator to maintain constant surface 
temperature. The primary limitation of the approaches in [1-5] is that the device employed is 
bulky and cannot be readily integrated with a microdevice. Few research groups have tried to 
control the temperature gradients at microscale by using active control strategies. For example, 
Hsieh et al. [8] and Wang et al. [9] employed an array type of microheater with an active 
compensation unit. This is in contrast with the conventional block heater employed to achieve 
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constant wall temperature for PCR application. Persat and Santiago [10] employed an 
unconventional way of chemical cycling (reducing the melting temperature of DNA with 
addition of solvent) which reduces thermal cycling. Wu et al. [11] fabricated a microheater and a 
thermal sensor inside a PDMS microchannel using injection molding and demonstrated 
uniformity in temperature distribution with the help of high precision control system. Although 
they succeeded in their objective, the level of complexity associated in their approaches cannot 
be ignored. This defeats the primary reason of employing a microdevice in the first place. In this 
paper, we propose a much simpler solution to this problem. 



2. PROPOSED DESIGN 

An effect that becomes particularly relevant at the microscale is conduction in the walls (or 
conjugate effect), which leads to a redistribution of heat flux and temperature in both the flowing 
fluid and the solid walls [12-15]. In most of the conventional flow situations, the cross section 
through which fluid flows is much larger as compared to the wall thickness. The primary mode 
of heat transfer is convection and the magnitude of axial conduction within the solid is small as 
compared to convection through the fluid. However at the microscales, the characteristic 
dimensions of the microchannel are of the same order or may be smaller than the thickness of the 
substrate. When such a situation arises, the axial conduction through the solid cannot be 
neglected; heat transfer then happens by convection in the fluid coupled with axial conduction in 
the solid. 

In several applications the bottom wall of the substrate is subjected to a constant heat flux 
boundary condition. If the conjugate effects come into the picture then heat flux at the bottom of 
the microchannel is no longer constant (Fig. la). This redistribution in heat flux alters the surface 
temperature, resulting in a varying wall to fluid temperature difference (AT); as against a 
constant AT for uniform heat flux boundary condition (Figs. lc,d). It is noteworthy that heat flux 
can be further redistributed by employing geometry of varying cross section (Fig. lb). The idea 
presented in this work makes use of these two effects (conjugate and varying cross section) in 
order to reduce the surface temperature gradient in a diverging microchannel. 
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In the inlet part of the diverging microchannel the heat flux is high which leads to a higher 
temperature in this region (Figs. lc,d). As the fluid moves downstream, the heat flux decreases 
steadily and the rate of heating also goes down. This results in relatively lower temperature rise 
in the downstream section. The conjugate effects lead to an axial back conduction which tends to 
reduce the temperature gradients further. So, in effect, the temperature difference between the 
inlet and outlet is reduced to a great extent leading to a nearly isothermal boundary condition in a 
diverging microchannel. It is worthwhile to underscore that the proposed process is particularly 
advantageous in the sense that it utilizes a constant heat flux boundary condition which is 
practically easier to maintain but results in a near isothermal surface condition. 

To prove this concept, heating in diverging microchannel is undertaken with DI water as the 
working fluid. The bottom of the substrate is subjected to a constant heat flux boundary 
condition. A comprehensive study of the ensuing temperature distribution at the substrate bottom 
of the diverging microchannel is carried out experimentally as well as numerically (Section 4). A 
detailed parametric study has also been performed at various geometrical, thermo-physical and 
input conditions in order to understand their effect on the bottom wall boundary condition 
(Section 5). A model is also proposed to obtain the value of the surface temperature as a function 
of various governing parameters (Section 6). Finally, the design is demonstrated for PCR 
application (Section 7). 



3. METHODOLOGY EMPLOYED 

A brief description of the experimental and numerical methodologies employed in the present 
work is presented in this section. 

3.1 Experimental Investigation 

Figure 2 presents a schematic of the experimental setup employed in the measurements [16-18]. 
The experimental setup consists of separate water and electrical circuits. The pump used for 
pushing DI water through the microchannel provided volumetric flow rate ranging from 0.1 to 6 
ml/min. A pressure gauge equipped with data acquisition system recorded the pressure drop 
across the microchannel. A similar data logging arrangement was employed for recording 
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temperatures. The temperatures were measured by thermocouples placed at the inlet and outlet of 
the microchannel as well as on the bottom surface of the microheater. 

The microchannels employed for the purpose of experimental investigation were fabricated on 
silicon wafer using MEMS techniques, involving photo lithography and wet etching; using in- 
house fabrication facilities. The detailed fabrication process is discussed elsewhere [17]. The 
values of different geometrical parameters of the microchannel along with the associated 
uncertainties are provided in Table 1 . 

3.2 Numerical Investigation 

In order to obtain additional local information of the relevant parameters, three-dimensional 
numerical simulations were performed. For this purpose a microchannel having the same 
dimensions as that of the experimental test section was modeled numerically. A straight section 
having length of 5 mm was placed before the inlet to ensure that the flow is hydrodynamically 
fully developed before entering the microchannel. The silicon substrate thickness and length of 
the microchannel were taken as 286 |im and 20 mm respectively. The 1 |im thick microheater 
was also modeled. Figure 3 shows the microheater-substrate-microchannel assembly; the 5 mm 
straight reservoir towards the bigger width size is however not shown in the figure. 

The numerical simulations were performed using commercial package Fluent. The variation of 
fluid viscosity and thermal conductivity (k) with temperature (T) have been considered in the 
simulations. These are modeled as per the equations given below [19]: 

jU = cT d for300K<T<363 K (1) 

where c = 1.1 x 10 11 and d = -5.7. Note that the viscosity of water decreases by about 69% with 
increase in temperature, over the temperature range considered. The thermal conductivity 
however increases linearly with increase in temperature, and was modeled as 

k = a + bT for300K<T<373 K (2) 

where a = 0.0033, b = 0.0019. The increase in thermal conductivity is 11% over the temperature 
range considered. 
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Semi Implicit Method for Pressure Linked Equations (SIMPLE) algorithm was employed to 
solve the momentum and energy equations. Second order discretization scheme was employed 
for pressure equation; the momentum and energy equations were discretized using Quadratic 
Interpolation for Convection Kinetics (QUICK). A grid independence test was systematically 
carried out and accordingly 450,000 cells were selected for the simulations. The Nusselt numbers 
from the simulations were validated against the experimentally obtained value for a number of 
cases; the two agreed within 7%. Further validation of local temperature values is discussed 
through Figs. 5 and 6 presented later. The simulations are therefore deemed to be validated. 



4. PROOF OF THE PROPOSED DESIGN 

The combined results from experimental and numerical investigations will be presented and 
discussed in this section. First, the surface temperature distribution in an 8° diverging 
microchannel was compared with a straight microchannel having a uniform width of 1133 |um, 
depth of 86 |im, same solid thickness and material. Note that the width of uniform cross section 
microchannel has been selected based on the concept of equivalent hydraulic diameter [16]. The 
simulation was carried out at a mass flow rate of 8.33 x 10" 5 kg/s (5 ml/min) with a constant heat 
flux of 4.8 W/cm 2 supplied at the bottom wall. Figure 4 shows that the surface temperature 
difference between the inlet and outlet was significantly higher for straight microchannel (11 °C) 
as compared to a diverging microchannel (1 °C). In the case of diverging microchannel, the heat 
flux magnitude at the inlet is higher since the surface area at the inlet is smaller (Fig. 4b). This 
leads to an increase in the surface temperature near the inlet of diverging microchannel. The 
surface area of the microchannel keeps on increasing as we move downstream which leads to a 
decrease in the heat flux absorbed. Hence, the surface temperature near the outlet side is reduced 
to some extent. Also, the temperature gradients in the solid created by the flowing liquid lead to 
an axial back heat conduction which tries to reduce the temperature gradient. We therefore find 
that the temperature rise in a diverging microchannel is much less than a straight microchannel 
under identical conditions. This result is further verified through experiments. 

In order to make experimental measurements, the surface temperature was measured using a 
thermocouple positioned at different axial locations; and the temperatures were recorded under 
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steady state condition. The process was repeated over the entire length of the microchannel. 
Figure 5 shows the axial temperature variation for four different flow rates between 4.17x 10" 5 - 
9.17x 10" 5 kg/s. The reduction in surface temperature gradient is particularly evident at higher 
flow rates. Therefore, for the present microchannel configuration it is possible to reduce the 
surface temperature variation below 2° C over a large part of the microchannel by maintaining 
9.17x 10" 5 kg/s of water flow through a diverging microchannel. 

In addition to the axial variation, lateral variation of the temperature was also measured at both 
inlet and outlet locations. From Fig. 6, it can also be seen that the surface temperature variations 
along the lateral direction were not significant; and the temperature can be more or less assumed 
to be constant along the transverse direction. The reason for such uniformity of temperature in 
the transverse direction is the fairly high thermal conductivity of the silicon substrate (148 W/m- 
K). Figures 5 and 6 also present comparison of experimental and numerical results; the two are 
clearly in good agreement. 

The above results point to the fact that the surface temperature gradients reduce significantly for 
diverging microchannels. It can therefore be concluded that by judiciously varying the 
thermophysical properties and input parameters, a nearly isothermal wall boundary condition can 
be obtained for a diverging microchannel. In the next section, a detailed parametric study 
analyzing the effect of different geometrical, flow and thermo-physical properties on the wall 
boundary condition is undertaken to establish the above finding over a larger parameter range. 



5. EFFECT OF OTHER PARAMETERS 

A number of parameters affect the magnitude of axial back conduction and alteration in the heat 
flux distribution in the microchannel. These parameters are: geometrical (angle, length, depth, 
solid-to-fluid thickness ratio), thermo-physical (thermal conductivity ratio of solid-fluid), and 
input conditions (mass flow rate, heat flux). A detailed parametric study was carried out, as 
discussed under the subsequent subsections. 
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5.1 Effect of angle 

The three angles considered here are 1, 4 and 8°. The inlet width was fixed as 200 |iim and 
simulations were performed for a mass flow rate of 5x 10" 5 kg/s with bottom of the substrate 
maintained at a constant heat flux of 4.8 W/cm . The idea of redistribution in flux and its effect 
on obtaining isothermal conditions has been demonstrated successfully in Section 4 for 8° 
diverging microchannel. As the angle was increased from 1° onwards, there was relatively 
sharper increase in the area. So, the heat flux absorbed along the length of the microchannel as 
well as the Reynolds number gradients (thus Nusselt number) decrease steeply. Therefore, with 
increase in angle, the surface temperature gradients reduced to a great extent as evident from Fig. 
7. For example, surface temperature gradients (AT S /L; where (AT S is the difference in surface 
temperature and L is the microchannel length) for 1, 4 and 8° microchannel were observed to be 
0.6, 0.35 and 0.3 °C/mm respectively. 

The results therefore suggest that an increase in the angle leads to a reduction in the surface 
temperature gradient, and is therefore preferable. However, there is an upper limit to the angle 
(16°), beyond which flow reversal will come into the picture [16]; flow reversal would 
significantly alter the flow in the microchannel and the associated temperature and heat flux 
distributions. 

5.2 Effect of depth 

Effect of microchannel depth was also analyzed by simulating for three different microchannel 
depths - 86, 120 and 200 |im; while keeping the solid to fluid thickness ratio constant along with 
rest of the geometrical parameters. Again the flow rate was maintained at 5x 10" 5 kg/s and the 
wall heat flux as 4.8 W/cm 2 . Increasing the depth and maintaining a constant solid to fluid 
thickness ratio leads to a decreased Reynolds number (Re), which reduces the convection 
strength. As a result, the average surface temperature increases with increase in the microchannel 
depth, as can be observed from Fig. 8. Also, the surface temperature gradient is observed to 
decrease with depth. It can be observed that the gradients are 0.3, 0.2 and 0.1 °C/mm respectively 
for 86, 120 and 200 |im microchannels. This is because the reduction in Re with depth is more 
towards the narrower region. Therefore, the convection effect drops largely towards the narrower 
end, which leads to an increase in the surface temperature. On the contrary, towards the larger 
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end the drop in Re is not much and hence an increase in the surface temperature at this end is 
small which brings in the constant wall temperature condition. 

The microchannel depth therefore has an important effect on the temperature distribution; 
choosing a larger depth is advantageous from the point of view of maintaining a constant surface 
temperature. 

5.3 Effect of length 

In the next set of simulations, the effect of microchannel length on the ensuing surface 
temperature gradient was analyzed. The microchannel lengths considered were 10, 20 and 30 
mm with all other geometrical parameters remaining the same. The mass flow rate was 
maintained at 5x 10" 5 kg/s and the bottom wall was supplied with a constant heat flux of 4.8 
W/cm 2 . Figure 9 shows a nearly isothermal surface condition (-1 °C) for 10 mm length 
microchannel. It was further observed that as the length reduced, the surface temperature 
gradient also reduced. For 10, 20 and 30 mm length microchannels, the temperature gradients 
were respectively 0.10, 0.30 and 0.53 °C/mm. The reason behind this nature is the decrease in the 
convection effect with length. The decrease in Reynolds number for 30 mm length microchannel 
is substantially higher than that of 10 mm length microchannel. As a result the 30 mm length 
microchannel shows a larger surface temperature gradient owing to reduced strength of 
convection in the downstream direction. 

This result leads to a significant conclusion that a shorter length leads to a more uniform 
temperature distribution. Having a uniform temperature over a small length is particularly useful 
for microde vices. 

5.4 Effect of solid to fluid thickness ratio 

Conduction in solid becomes significant when the solid to fluid thickness ratio exceeds unity 
[12-15]. In order to examine the effect of solid to fluid thickness ratio on the surface temperature, 
the base case having t s /tf = 2.32 was compared against t s /tf = 1.5 and t s /tf = 4. MicroChannel depth 
was maintained constant as 86 |im while the solid thickness was varied appropriately to obtain 
the desired ratio. This was done while maintaining all other geometrical parameters same as the 
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base case. Again, the heat flux was maintained at 4.8 W/cm 2 and the mass flow rate was kept 
5 x 10" 5 kg/s for this set of simulations. 

It was observed that the surface temperature gradients reduce when t s /tf increases owing to an 
increase in the strength of axial back conduction due to the increase in solid cross sectional area 
(Fig. 10). However for the range of values considered here, the difference in surface temperature 
gradient is insignificant. A drastic variation of t s /tf is required to observe substantial change in 
the surface temperature distribution. 

5.5 Effect of solid to fluid thermal conductivity ratio 

Simulations were also performed to study the effect of solid conductivity (k s ) to fluid 
conductivity (k£) ratio with working fluid fixed as water. For this purpose, the base case was 
taken as silicon (k s /kf = 247) and the results compared against steel (k s /kf = 27) and copper (k s /kf 
= 646) as the microchannel substrate materials. The heat flux was again maintained as 4.8 W/cm 2 
and the mass flow rate as 5x 10" 5 kg/s. Since conductivity ratio is the only variable, it affects the 
intensity of axial back conduction. Increasing the conductivity ratio reduces the thermal 
resistance (L/k s A; where L is the microchannel length and A is the cross section area) to axial 
back conduction. Therefore, the average surface temperature reduces for silicon and copper as 
compared to steel (Fig. 11). The temperature gradients come out to be 0.35, 0.3 and 0.28 °C/mm 
for steel, silicon and copper, respectively. 

The above results suggest that unless the conductivity ratio is varied drastically, a substantial 
effect on the surface temperature gradients will not be seen. That is, the results presented in this 
study using silicon wafer is equally applicable to many other solid-liquid combinations as well. 
Any material which does not react with the biological sample and has a good thermal 
conductivity can be employed to make microchannels, making this methodology all the more 
attractive. 

5.6 Effect of mass flow rate 

Additional simulations were performed for 8° diverging microchannels in order to observe the 

effect of mass flow rate on the temperature distribution at the bottom surface. The mass flow 

rates considered were 2, 3 and 5 ml/min at supplied heat flux of 4.8 W/cm keeping all other 
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geometrical parameters same. It is observed that the surface temperature gradient as well as the 
average surface temperature reduced substantially at higher flow rates (Fig. 12). The temperature 
gradients for 2, 3 and 5 ml/min are 0.45, 0.3 and 0.05°C/mm respectively. At higher flow rates, 
the convection heat transfer becomes substantially higher than the wall conduction. Therefore the 
relative magnitude of conjugate effect is reduced with an increase in the flow rate. 

Hence, both lowering of the average surface temperature as well as the reduction of surface 
temperature is brought about at higher flow rates. However, employing higher flow rate would 
lead to a higher pressure drop penalty and should therefore be avoided. The constancy of the wall 
temperature should rather be maintained by suitably altering some other parameter. 

5.7 Effect of heat flux 

Simulations were also performed in order to study the effect of heat flux supplied at the bottom 
wall. The mass flow rate of 5xl0" 5 kg/s was passed through a 8° diverging microchannel at three 
different heat fluxes of 2.4, 4.8, 9.6 W/cm 2 . The average temperatures went down with reduction 
in heat flux as expected (Fig. 13). Along with that, the temperature gradients were observed to be 
0.1, 0.3 and 0.65 °C/mm for 2.4, 4.8 and 9.6 W/cm 2 , respectively. Again a near constant surface 
temperature was observed for lower heat flux. That is, a higher heat flux leads to a higher 
temperature in addition to a higher temperature gradient. 

Combining the results of this section with Section 5.3, it is noted that a lower heat flux over 
smaller segments of length is better than a higher heat flux over a larger length from the point of 
view of maintaining a nearly isothermal surface condition. 

5.8 Summary of the parametric study 

The experimental and numerical results showed that the combination of varying cross-section 
and conjugate effects lead to a redistribution of heat flux which significantly alters the wall heat 
flux boundary condition. This also affects the temperature distribution in the microchannel. A 
detailed parametric study was carried out to study the effect of different geometrical, thermo- 
physical and input parameters on the wall boundary condition which also provides an idea about 
the relative strength of convection and wall conduction in diverging microchannel. 
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Table 2 presents a summary of the parametric study. It is observed that the temperature gradient 
(AT S /L) decreases asymptotically with divergence angle (Case I). This could be attributed to a 
reduction in Reynolds number and heat flux at inner walls in the downstream direction with 
increase in angle. In Case II, it is observed that the average surface temperature increases with 
depth of microchannel owing to a decrease in Re. However, AT S /L is observed to reduce with 
depth since the change in Re is more towards the narrower region as compared to the wider 
region. In Case III, the length of the microchannel was varied and it was found that a shorter 
length microchannels exhibits better surface temperature homogeneity (lower AT S /L). The reason 
for this is larger convection in shorter length microchannel (Re is higher) than in longer one with 
same angle of divergence. Also, there is a reduction in the resistance to the axial back conduction 
with decrease in length of microchannel. In the next two cases (Cases IV and V), the effect of 
conjugate heat transfer on AT S /L was studied. In this, solid to fluid conductivity ratio and 
thickness ratio were varied and it was concluded that only a drastic variation in these parameters 
would have a significant bearing on the surface temperature variation. Furthermore, the effects of 
input parameters were studied. It was found that increasing the mass flow rate reduced the 
temperature gradient (AT S /L). After this, simulations were performed by varying the substrate 
bottom wall heat flux and greater temperature homogeneity was observed when the heat fluxes 
were lower in magnitude. 

In view of the detailed parametric study, the design parameters to achieve a nearly isothermal 
condition can be selected within the limit prescribed for the important parameters as: divergence 
angle (8-16°), depth (> 86 \±m) and length (< 20 mm) of microchannel along with mass flow rate 
(>3 ml/min) and heat flux (< 4.8 W/cm ). Hence, it can be concluded from the parametric study 
that under the set of optimal parameters, diverging microchannel would yield a near isothermal 
boundary condition (±0.5 °C). 



6. MODEL FOR ESTIMATING THE SURFACE TEMPERATURE 

This section presents a model for estimating the surface temperature as a function of the various 
governing parameters. The role of various governing parameters in the direction of achieving 
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isothermal wall condition has been identified in the previous section. Accordingly, a diverging 
microchannel with divergence angle: 0, depth: H and length: L is considered. 

Consider a control volume of diverging microchannel as shown in Fig lb. Using energy balance 
over the control volume we can write at steady state: 

E =E , (3) 

in out V / 

where, E in and E out are respectively the net energy entering and leaving the control volume. 

Now, based on the above parametric study, it is assumed that most of the energy gets convected 
with the flow. Therefore, we can write Eq. 3 as, 

dq (x)P(x)dx = mc p dT f , (4) 

where, q (x)is heat flux supplied (W/m 2 ), P(x) is the perimeter of diverging microchannel (m), 
m is the mass flow rate (kg/s), c p is the specific heat capacity (J/kg-K), and T f is local fluid 
temperature (K). Integration of Eq. 4 can be shown to give 

T f (x) = T fi +^—x (5) 
mc p L 

where and q are the inlet temperature (K) and power (W) supplied to the diverging 
microchannel of length L (m). Eq. 5 gives the variation of local fluid temperature T f (x) . 

Now consider a nearly isothermal wall boundary condition, Eq. 4 can be re-written as 

dT f 

rhc p —t = h(x)P(x)AT , (6) 
dx 

where h(x) is the convective heat transfer coefficient (W/m 2 -K) and Aris the difference 
between surface and fluid temperatures, i.e. (T s - T f ^in K. For nearly isothermal wall boundary 

d(AT) dT f 

condition we assume that AT varies with the same rate as 7}, i.e., -« — -. (This 

dx dx 

assumption is evaluated towards the end of the section.) Therefore, Eq. 6 becomes; 
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d(AT)_ h(x)P(x) 



AT 



dx , 



mc„ 



(7) 



Eq. 7 can be integrated from inlet to outlet assuming that h is a constant calculated at the 
characteristic location for diverging microchannel [16]: 



7d(AT) h f 
[ '- = — -\P(x)dx 



AT; 



mc 



(8) 



p o 



The heat transfer coefficient can be obtained from a known correlation for Nusselt number: 



Nuk f aRe^ k f 
h= f - f 



(9) 



where a = 0.301 and /? = 0.563 are obtained from the curve fitting of experimentally derived 
Nusselt number versus Re curve (not shown). In the above equation, kf is the thermal 
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conductivity of fluid (in W/m-K) and Reynolds number can be calculated from Re = — - . The 

A/j 

4A 

hydraulic diameter D h = — where, A, P are the cross sectional area and perimeter of diverging 

microchannel calculated at the characteristic location. The characteristic location lies at L/3 from 
the inlet [16]. 

After some algebra, Eq. 8 yields 



{Ts, 0 -T SJ ) = {T SJ -T fi )ex V 



hL 



mc 



w. +Ltan 



p v 



+ 2H 



T —T - 

1 f,i 



mc 



p J 



(10) 



from which the average surface temperature can be obtained as 



s 2 



«Re^ k f PLj 



1 

+ — 
2 



exp 
qD l 



aRe fi k f L 



C 



D h™ c P 



W; + L tan 



v 



+ 2H 

J 



2T fi + V + ■ 

} ' aRe fi k f PL mc p 



(11) 
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in terms of the relevant governing parameters. 

The surface temperature from Eq. 11 has been validated against experimental data in Table 3. 
The model clearly predicts the experimental data in the various cases tested to within 1 °C. The 
slightly higher value predicted by the model could be because we have neglected the effect of 
conduction heat flux in the model. 

The above model, which has a strong theoretical basis, is therefore reasonably accurate. There 
are two main assumptions involved in the derivation of the model: (i) AT varies with the same 
rate as 7}, and (ii) h calculated at the characteristic location is representative over the entire 
microchannel. These assumptions are required to avoid unnecessary mathematical complexity 
without comprising the accuracy of the predictions. The effect of these assumptions is apparently 
small, at least for the parameter range considered here, as evident from the comparison in Table 
3. We believe that these assumptions are not too severe even otherwise. 

The model can be used to predict the temperature that the surface will attain for a given set of 
parameters. The model can further help in deciding the optimal design parameters of diverging 
microchannel in order to achieve a nearly isothermal boundary condition. The model was indeed 
utilized in this latter manner for PCR application as presented in the next section. 

7. APPLICATIONS AND SIGNIFICANCE OF THE DESIGN 

In certain biological applications where maintaining a constant temperature (typically 37°C) is of 
utmost importance in order to keep the cells alive, there is a need to design portable systems 
which can maintain the desired temperature [1, 2]. Our results show that this is possible and the 
desired temperature can be maintained within ±0.5 °C. In order to achieve the same objective, 
researchers [6,7,11] have proposed ingenious albeit intricate ways. The appeal of the proposed 
method lies in its simplicity. By allowing a liquid to flow through a diverging microchannel 
supplied with constant heat flux at its bottom wall, a constant temperature can be achieved at the 
microchannel walls. The proposed design therefore adds portability to the existing processes by 
eliminating the need of an incubator. This also suggests that for a given geometry and properties 
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of diverging microchannel, any constant wall temperature can be obtained by changing the input 
parameters appropriately (i.e. with the help of Eq. 11). 

We now demonstrate our concept on a more complex example of PCR. The objective here is to 
maintain the different stations at the required temperatures of 95, 55, 72 °C (as explained in 
Section 1) so that PCR can be effectively carried out on a single chip. Three thermally insulated 
differently optimized diverging microchannels connected in series are employed for this purpose 
(Fig. 14a). The microchannels have different depths and different input heat fluxes. All the 
geometrical and input parameters are selected based on our model (Eq. 11). The results are 
obtained numerically and presented in Fig. 14b. The results are extremely encouraging in that 
three constant temperatures (-95, 57 and 73 °C) can be achieved in the three diverging 
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microchannels with heat flux of 10.5, 0.09 and 3.8 W/cm , for mass flow rate of 5xl0" J kg/s. 
The proposed design is therefore an important step forward in successful implementation of 
various lab-on-chip technologies. 



8. CONCLUSIONS 

In this study, we show that microchannel of varying cross section brings about redistribution in 
heat flux; also wall conduction (or conjugate effect) comes into picture at microscales. These 
effects can be exploited to our advantage to obtain a constant wall temperature surface condition. 
Heating in microchannels was studied numerically and experimentally in order to observe the 
effect of redistribution in heat flux and the ensuing surface temperature gradient. A comparative 
study performed between straight and 8° diverging microchannels showed a significant reduction 
in the surface temperature gradient for the latter case, lending credence to our proposition. 

The proposed idea was also verified experimentally in an 8° diverging microchannel over a wide 
range of flow rates. The lateral distribution of temperature showed an insignificant variation 
while the axial temperature gradient reduces substantially; thus implying that the surface 
temperature is almost constant. A detailed parametric study of the effect of various geometrical 
(angle, depth, length and solid-to-fluid thickness ratio), thermo-physical (thermal conductivity of 
solid) and input parameters (mass flow rate and heat flux) on the microchannel surface was 
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undertaken to explore this finding over a wider parameter range. The outcome of this parametric 
study suggested the selection of important parameters as: divergence angle (8-16°), depth > 86 
|Lim, length < 20 mm; along with mass flow rate > 3 ml/min and heat flux < 4.8 W/cm . For this 
set of optimal parameters, diverging microchannel would yield a nearly isothermal boundary 
condition with temperature difference between inlet and outlet of less than 1°C. 

A mathematical model is also proposed to obtain the optimal geometrical parameters for a given 
set of input parameters (mass flow rate and heat flux). Equation 11 derived in the paper can be 
employed for the purpose of designing the microchannel. 

This constant temperature boundary condition in a diverging microchannel can be effectively 
utilized in certain biological applications where the samples are highly sensitive to temperature 
fluctuations. Finally, a successful demonstration of the proposed design is shown numerically for 
an on-chip PCR. The simple and novel approach presented here can be employed in several other 
applications as well. 
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Table 1 : Geometrical specification of the microchannel employed in the present study 



Feature 


Size 


Uncertainty 


Units 


Smaller Widtn (Wi) 


Z/U 


+ z 


|im 


i^arger wiQtn (,w 0 j 


JUJZ 


-t- z 


|im 


Length (L) 


20 


-0.1 


mm 


Height (H) 


86 


+ 1 


|im 


Divergence Angle (0) 


8 


0.5 


0 


Hydraulic Diameter (Dh) 


157 


1.33 


|im 
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Table 2: Summary of Parametric study showing dependency of AT/L on each parameter. 



(Case) 
Parameters 


Wall 
Temperature 
gradient 


Non dimensional Geometrical 
Parameters 


Input Parameters 


Remarks 
























AT S /L 
(°C/mm) 


9 


w/H 


L/H 


ts/H 


ks/kf 


m 

(ml/min) 


q" 

(kW/m 2 ) 








0.875 


0 


2 


233 


2.33 


247 


3 


48 






(I) Effect of 


0.600 


1 


4 


233 


2.33 


247 


3 


48 


Strong 




Angle 


0.350 


4 


8 


233 


2.33 


247 


3 


48 


function 






0.300 


8 


13 


233 


2.33 


247 


3 


48 






(II) Effect of 
depth 


0.300 
0.200 


8 
8 


13 
9 


233 
167 


2.33 
2.33 


247 
247 


3 
3 


48 
48 


Strong 
function 


"onvection 




0.100 


8 


6 


100 


2.33 


247 


3 


48 






(III) Effect 
of length 


0.100 
0.300 


8 
8 


8 

13 


116 

233 


2.33 
2.33 


247 
247 


3 
3 


48 
48 


Strong 
function 






0.533 


8 


19 


349 


2.33 


247 


3 


48 






(IV) Effect 
of thickness 


0.300 
0.300 


8 
8 


13 
13 


233 
233 


1.50 
2.33 


247 
247 


3 
3 


48 
48 


weak 
function 


o 




0.300 


8 


13 


233 


4.00 


247 


3 


48 




© 


(V) Effect 
of Thermal 
conductivity 


0.275 
0.300 
0.350 


8 
8 
8 


13 
13 
13 


233 
233 
233 


2.33 
2.33 
2.33 


646 
247 
27 


3 
3 
3 


48 
48 
48 


Moderate 
function 


= 

CTQ 
P 
p+ 
fD 


(VI) Effect 
of Mass 
flow rate 


U.4j(J 
0.300 
0.050 


8 
8 
8 


13 
13 
13 


233 
233 
233 


2.33 
2.33 
2.33 


247 
247 
247 


2 
3 
5 


48 
48 
48 


Strong 
function 


o 

o 

< 


(VII) Effect 
of Heat flux 


0.075 
0.300 


8 
8 


13 
13 


233 
233 


2.33 
2.33 


247 
247 


3 
3 


24 
48 


Strong 
function 


CD 
O 
r+ 

S" 




0.650 


8 


13 


233 


2.33 


247 


3 


96 







20 



Table 3: Validation of proposed model with experimental data (T s denotes average surface 
temperature). 



Mass flow 
rate (kg/s) 


Power 
Supplied 
\ vv ) 


T ,°C 
(Experimental) 


T s , "C (Proposed 
model Eq. 11) 


\T s \ -t\ ), 

\ 1 model \expt J 

°C 


3.33 x 10" 5 


3.5 


49 


50 


1 


5.00 x 10" 5 


3.5 


44 


45 


1 


6.67 x 10" 5 


3.5 


41 


42 


1 


8.33 x 10"' 


3.5 


39 


41 


2 


10.0 x 10" 5 


4.0 


40 


41 


1 




q TT a bs at inner wall with 
CHT pi verging) (Fig. lb) 

qf' at inner wall with 
CHT (Fig. la) 



<l" abt at inkier wall 
vathout CHT 



\ 0". ™ at outer v.- all 



A 

I 



X 

(c) 

T s uc = Surface temperature in uniform cross section 
T s .dc — Surface temperature difference in diverging 



.^t - *■ 



- " T 



s.uc 



X 
(d) 

Figure 1: Schematic representation of the proposed method bringing about constant wall 
temperature condition, (a) and (b) show the heat flux at the walls for uniform cross section and 
diverging microchannel respectively. A graphical representation of the heat flux variation and 
temperature variation along the flow direction are provided in (c) and (d) respectively. (CHT 
stands for conjugate heat transfer.) 
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Figure 2: Schematic of the experimental setup 
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Figure 3: Three-dimensional model employed in the numerical simulations. A diverging (8°) 
microchannel with t s /tf = 2.32 (not to scale) is shown here. 
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Figure 4: Comparison of (a) surface temperature and (b) heat flux distribution in a diverging and 
a straight microchannel. 
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Figure 5: Surface temprature variation along the flow direction in a diverging microchannel at 
q" = 4 W. (Symbols represent experimental data whereas lines show the corresponding 
numerical result.) 
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Figure 6: Comparison of lateral surface temperature variation at (a) inlet and (b) outlet in a 
diverging microchannel, for m = 4.17xl0" 5 kg/s and q sup = 4 W. (Symbols represent 
experimental data whereas lines show the corresponding numerical result.) 
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Figure 7: Effect of angle on surface temperature variation in diverging microchannel. 
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Figure 8: Effect of depth variation on surface temperature distribution in diverging microchannel 
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Figure 9: Effect of length on surface temperature distribution in diverging microchannel. 
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Figure 10: Effect of t s /tf variation on surface temperature distribution in diverging microchannel. 
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Figure 1 1 : Effect of solid to fluid conductivity ratio on surface temperature variation in diverging 
microchannel. 
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Figure 12: Effect of mass flow rate on surface temperature distribution in diverging 
microchannel. 
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Figure 13: Effect of heat flux variation on surface temperature distribution in diverging 
microchannel. 
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Figure 14: (a) Schematic representation of a microdevice for on-chip PCR process. The 
microdevice consists of three diverging microchannels of different geometrical configurations, 
(b) The temperature profile in the three stations (diverging microchannels). 
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